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Abstract: It has recently been pointed out that simple scaling properties of Polyakov 
correlation functions of gauge systems in the confining phase suggest that the ratios of k- 
string tensions in the low temperature region is constant up to terms of order T^. Here we 
argue that, at least in a three-dimensional Z4 gauge model, the above ratios are constant 
in the whole confining phase. This result is obtained by combining numerical experiments 
with known exact results on the mass spectrum of an integrable two-dimensional spin model 
describing the infrared behaviour of the gauge system near the deconfining transition. 
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1. Introduction 

In most confining gauge theories, besides the fundamental string (of tension cr) which is 
formed between a pair of static sources in the fundamental representation /, there is the 
freedom of taking the sources in any representation Si. 

If, for instance, the gauge group is SU(iV) there are infinitely many irreducible rep- 
resentations at our disposal. However, as the sources are pulled apart, no matter what 
representation is chosen, the asymptotically stable string tension depends only on the 
A^— ality k of i^, i.e. on the number (modulo A^) of copies of the fundamental representa- 
tion needed to build Si by tensor product, because all representations with the same k can 
be transformed into each other by the emission of a proper number of soft gluons. As a 
consequence the heavier strings decay into the string of smallest string tension cr^. . The cor- 
responding string is referred to as a k-string. This kind of confining object can be defined 
whenever the gauge group admits more than one non trivial irreducible representation. 

Much work has been done in the study of k-string tensions in the continuum pl-pl as 



well as on the lattice l0-[16l 
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In a previous work [0] , some of us have argued from simple scaling properties of suit- 



able Polyakov loop correlators that these string tensions have the following low temperature 
asymptotic expansion 



where c is the central charge of the underlying 2D conformal field theory describing the IR 
behaviour of the k-string. As a consequence, their ratios are expected to be constant up 
to terms: 



The low temperature data presented in support of this expectation were taken from 
Monte Carlo simulations on a particular system, namely a (2+1) -dimensional Z4 gauge 
model, which is the simplest exhibiting more than just the fundamental string. 

The main conjecture we want to verify in this work is that ak{T)/a{T), at least in 
that Z4 gauge system, is in fact independent of the temperature in the whole of the con- 
fining regime. To check this idea, a handy fact comes useful, namely that, as the system 
approaches the deconfinement transition, and the string picture begins fading, another 
approach is made available by the Svetitsky-Yaffe (SY) conjecture [|l8|, which allows to 
reformulate the system in a totally different perspective, based on a two-dimensional in- 
tegrable theory in which, however, the near-Tc counterpart of the low-temperature result 
cited above can be nicely found. 

It turns out that the deconfinement transition of the 3D Z4 gauge model is second 
order and, according to the SY conjecture, belongs to the same universality class of the 
2D symmetric Ashkin- Teller (AT) model. As a matter of fact, such a model possesses 
a whole line of critical points along which the critical exponents vary continuously. The 
SY conjecture tells us that if a (2+l)-dimensional gauge model with center Z4 displays a 
second-order transition, then its universality class is associated to a suitable point of the 
critical line of the 2D AT model. For instance, it has been argued [^9| that the critical 
(2+l)D SU{4:) gauge theory belongs to the universality class of a special point of the AT 
model, known as the four-state Potts model. More generally, the class of models with 
gauge group Z4 depends on two coupling constants a and /?, and the universality class of 
the deconfining point P varies with the ratio a/p. 

The two-dimensional AT model can be seen in the continuum limit as a bosonic con- 
formal field theory plus a massive perturbation (i. e. a Sine-Gordon theory) driving the 
system away from the critical line. Thus, a map between (a neighbourhood of) the AT 
critical line and the Sine-Gordon phase space is provided. 
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This theory is integrable, and the masses of its hghtest physical states (first sohton 
and first breather mode, of masses M and Mi) correspond to the tensions cr{T) and cr2(T) 
near Tc, whose ratio, in this context, can be analytically evaluated and turns out to be 

lim ^ = ^ = 2sin^(2.-l) , (1.3) 

T^Tc a{T) M 2^ ' ' ^ ' 

where is the thermal exponent in two dimensions. 

As a consequence, on the gauge side, we have two different ways to verify the conjecture. 
One is to directly estimate the ratio Mi /M by measuring the Polyakov-Polyakov correlators 
in the two non-trivial representations of Z4. The other is to evaluate the thermal exponent 
of the gauge system at the deconfining temperature. Either method gives a value of Mi/M 
which nicely agrees with the ratio (T2/1T evaluated at T = . 

1.1 The (2+l)D Z4 gauge model and its dual reformulation 

The most general form of Z4 lattice gauge model admits two independent coupling con- 
stants, with partition function 

in which the gauge field Ui on the links of a cubic lattice is valued among the fourth roots 
of the unity and the sum in the exponent is taken over the elementary plaquettes of the 
lattice. Such a theory can be reformulated as two coupled Z2 gauge systems: 

m,Pff) = ll E e^.[^H^^+^^)+/^//^^^^] ,{U, = llUi ; V, = l[VO. 

I {;7i=±i,y,=±i} lep lep 

(1.5) 

Prom a computational point of view it is useful to exploit a duality relation, switching 
to a spin model, that in this case is a 3D AT model, expressed as a double Ising spin field 
plus a coupling term between the Ising variables {a} and {r}: 

SAT{a,P) = - "^iPi^^xf^y + TxTy) + a{ax(JyTxTy)] ; (1.6) 

{xy) 



the duality is implemented by 

1 ^ |"(coth/3y + tanh/?j tanh/3jj)(coth /3j + tanh/3y coth. j3ff) 



2 + tanh Pff + coth Pff 



(1.7) 



4 



1 + tanh^ Pf tanh f3ff 



tanh (5f + tanh.Pff 



/?=-ln 2 ^ ^ . (1.8) 



The choice of working in the dual spin version of the gauge system is strongly motivated 
by the availability of highly efficient nonlocal Monte Carlo algorithms, in which, moreover, 
any kind of gauge-invariant observable can be directly embedded in the update procedure. 
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The Z4 gauge model admits, in addition to the fundamental string, a k-string with 
k = 2, corresponding to taking the sources in the double-fundamental representation f ® f ■ 

The phase diagram of this 3D model at T = has been studied long ago |2^, 21]. The 
deconfinement transition, which is weakly first order in the region we are interested in, 
becomes second order at finite T, therefore, according to SY conjecture, is described by 
the order-disorder transition of a 2D AT model. As anticipated in the Introduction, such 
a transition forms a whole 1-dimensional manifold of critical points (see Figure |^). Along 
this line, the critical indices (and other universal quantities as well) vary continuously, the 
endpoints representing a decoupled double Ising system and the 4-state Potts model. The 
choice of the point in the phase space in which to work is therefore a crucial issue. 



Ordered phase 



Symmetric phase 



J4 



Figure 1: Sketch of the phase diagram of the Ashkin- Teller model in two dimensions (see Eq. 2.1 
for the meaning of the constants J, J4). The point labelled / corresponds to a pair of decoupled 
critical Ising systems (J4 = 0, J = Jl^™^), while the P4 point (J4 = J = J^"'*") represents the four- 
state Potts critical point. The solid (self-dual) line between them is the one-dimensional critical 
manifold along which the critical indices vary with continuity. Beyond P4, the self-dual line is no 
longer critical, and another phase appears (bounded by the two dashed lines departing from P4). 
The point P indicates the image of the confined three-dimensional gauge system near according 
to the Svetitsky-YafFe conjecture. 



From the data in |17], obtained by means of finite-temperature measurements of 
Polyakov-Polyakov correlation functions, and particularly from those referring to the point 
P identified by (a,/3) = (0.050,0.207), the string tensions a and 02 can be evaluated in 
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the T ^ limit as temperature-independent quantities: 

ac? = 0.02085(10) , 

f72 = 0.03356(22) , (1.9) 

where a is the lattice spacing. Their ratio, which has been argued to equate the central 
charge of the CFT related to the 2-string, is then given by 

— = 1.610(13) . (1.10) 
a 

2. The Svetitsky-Yaffe conjecture and the Sine-Gordon model 

The mapping induced by the Svetitsky-Yaffe conjecture leads to a substantial simplification 
in the study of the critical properties of the deconfining transition, allowing to study it as 
a standard symmetry-breaking transition which takes place in a spin model. In the present 
case we deal with the symmetric Ashkin- Teller model in two-dimensions. 

The action for this model has the same form of Eq. ( |1.6D , but to clarify the fact that 
this is the 2-dimensional AT system reached via SY conjecture (as opposed to the three- 
dimensional one dual to the original gauge theory), we will relabel the spin fields and 
coupling constants with other names within this context: 

Sat = - + o-^o-y) + ■ (2-1) 

Such a model has been extensively studied in the past, and a number of analytic and 



numerical results have been discovered 22 1. Its phase diagram is exactly known (Fig. |l]). 



The critical line is characterised by the fact that it is self-dual and separates a disor- 
dered phase from an ordered one. Since the symmetry which is spontaneously broken by 
crossing it is the global Z4 symmetry of the model, this is also the critical manifold corre- 
sponding to the deconfining phase transition of the gauge model. Hence in the following 
we will concentrate on (a part of) this critical line. 

Important analytic results related to the critical line have been worked out by a direct 
solution of the lattice model |22|. Another interesting approach was carried on in the paper 



|23|] , by considering the c = 1 conformal field theory describing the critical line, i.e. the 
Gaussian model. In such a context it was possible to establish an exact correspondence 
between lattice/continuum operators and to compute exactly their conformal dimensions. 

The advantage of working in the field theoretical setting is that one has the possibility 
to study the off-critical behaviour near the critical line in a very natural way. This is of 
great convenience for us, since we will be ultimately interested in the mass spectrum of the 
theory in the high temperature phase. 
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Let us briefly sketch how to obtain a field theoretic description of the scaling region 
near the critical self-dual line of the model. Since the correlation length remains larger 
than the lattice spacing, the model can be described by the following QFT: 

At = Am^ + Am^ +r j (fx{ei{x)+e2ix)) +p J (f x ei{x) e2{x) , (2.2) 

where the meaning of such an expression is quite transparent when compared to the Hamil- 
tonian of the lattice model. The first two terms stand for the conformal field theories with 
c = 1/2 describing the critical behaviour of the two Ising models, and the latter two terms 
are respectively the relevant thermal perturbation (ej(x) are the energy operators in the 
two copies of the Ising model), and the marginal one which moves the system along the 
critical line. It is also evident by comparison that the couplings r and p are substantially 
a reformulation of the former J and J4 respectively. 

In other words, when r = 0, p ^ 0, the critical line is described by a compactified free 
massless boson (Gaussian model). A full analysis of the critical line of the Ashkin- Teller 



model by means of the Gaussian model was established in |23]. 

The general Hamiltonian with r 7^ describes the model outside the critical line. It is 
useful to note that it can be seen as a perturbation of the Gaussian model, and in such a 
bosonic language the thermal perturbation can be written as cos P(p, where /3 is a marginal 
parameter equivalent to p. Hence we are left with 

At = j <fx - T COS/?(/P^ , (2.3) 

which is the action of the Sine-Gordon model. In this notation, we have the high-temperature 
phase for r > 0, and the low-temperature one for r < 0. It was also shown that the previ- 
ous model describes the Ashkin- Teller model in the range 27r < /32 < 67r. Actually we are 
interested in the narrower range 27r < < Att fr om the critical 4-state Potts model to two 
decoupled critical Ising models. Furthermore, since the confined phase of the gauge theory 
is mapped in the high-T phase of the Ashkin- Teller model, we will only consider the case 
r > 0. 

Such a QFT is of particular interest because it is integrable, and this is the main reason 
for rewriting the action of the model near the critical point in a bosonic form^. 

Integrability means that an infinite number of integrals of motion exists. The main 
consequence in (l+l)-dimensions is the fact that the scattering theory is very constrained, 
because the S'-matrix is factorised in products of two-body interactions, and inelastic pro- 
cesses are forbidden. 



^It is always possible to fermionise the action of the Sine-Gordon model in order to obtain an integrable 
fermionic theory with the same scattering matrix, namely the massive Thirring model. 



-6- 



These facts allow to write down the so-called Yang-Baxter equations for the 2-particle 
S-matrix. Then, such an S-matrix can be computed exactly by imposing the previous equa- 
tions and the usual requirements of unitarity and crossing (for a review about Integrable 
QFTs see 11). 

An obvious consequence is that also the spectrum of the masses of the bound states of 
the theory is known exactly, since they are represented by the simple poles of the S'-matrix 
in the physical strip. 

It is worth to recall that an important consequence of the Svetitsky-Yaffe conjecture is 
that the ratio of a given string tension over the temperature of the gauge theory is mapped, 
near Tc, onto a corresponding mass of the spin model spectrum. Then, the ratio of string 
tensions 02! o is mapped onto a suitable ratio of masses of the Sine-Gordon model which 
is known exactly as a function of universal quantities. 

In the following we will make some quantitative considerations about the qualitative 
picture given above in order to make some predictions which will be useful in the context 
of the gauge theory. Since a detailed analysis of the properties of the scattering theory of 
the Sine-Gordon model in the context of the Ashkin- Teller model has been done in p5| , 
we will refer to those papers for the details. 

2.1 Operators correspondence, mass spectrum and correlation functions 

The last ingredient we need before exploiting the map to the Sine-Gordon model at its 
best is the correspondence between the Polyakov loops in higher representation and the 
operators of the Ashkin- Teller model. 

We already know from the Svetitsky-Yaffe original work that the Polyakov loop in the 
fundamental representation corresponds to the spin operator. Then, following the same 
reasoning used in |^], it is possible to deduce that the Polyakov loop in the double fun- 
damental representation is related to the so-called polarisation operator V = where 
and are the spin variables defined in ( |2.1| ). Its bosonic form and the corresponding 
anomalous dimensions are given by 

P = sin|v9, Xr = ^; (2.4) 

we also notice that {V) = in the high-T phase of the model. 

Sine-Gordon mass spectrum |]2^ : The exact knowledge of the S-matrix allows to access 
to the exact mass spectrum of the theory. Without entering into the details, the spectrum 
of the SG model is given by a soliton/anti-soliton doublet of fundamental particles of mass 
M, and a number of soliton/anti-soliton bound states, called breathers -B„, whose number 
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is a function of /3 . By defining the couphng constant ^ in the following way 



we have that for > vr, i.e. /J^ > An, no bound states are present and hence the spectrum 
is given by the soliton/anti-soliton doublet only (repulsive regime). 

For ^ < TT, i.e. /J^ < 47r, we are in the attractive regime and the breathers appear as 
simple poles of the S-matrix. Their number and masses are given by the following formula 



n 

M„ = 2M sin - ^ , 1 < n < 



(2.6) 



Since we are interested in the range 2tt < (3'^ < An, we immediately realise that, outside 
the two decoupled Ising point at /J^ = Air, we always have at least one breather of mass 
Ml (for 27r < < S/Stt we also have the breather M2 which is however irrelevant for our 
analysis) . 

The next step is to associate particle states to operators in the high temperature phase. 
It has been done in by taking into account their properties of symmetry and locality. 
The result is that the spin operator only couples to particle states with topological charge 
equal to one, i.e. involving an odd number of solitons (or antisolitons), and couples to 
even-labelled breathers only. On the contrary the polarisation operator couples to neutral 
particle states only, i.e. states with the same number of solitons/antisolitons and odd- 
labelled breathers. 

As a consequence the spin operator is naturally associated to the mass of the soliton, 
and the polarisation operator is associated to the mass Mi of the breather Bi . This means 
that the string tension of the Polyakov loop in the fundamental representation corresponds 
to the mass of the soliton, and the the string tension of the double fundamental corresponds 
to the first breather. Hence, following the Svetitsky-Yaffe conjecture, the ratio of string 
tensions in the confining phase near the transition is given by 

— — = 2 sm - . (2.7) 
M 2 ^ ^ 

This result, being a dimensionless ratio, is expected to be universal in the limit r — > 0. 
This fact can be explicitly seen by expressing the coupling ^ in terms of some critical 
exponent. Not surprisingly, it is indeed possible because the theory is solved also at the 
critical point in terms of the Gaussian model. By comparing the power-like behaviour of 
the energy-energy correlator of the Gaussian model: 

(e(x)e(O)) 

Gaussian 
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with that expected from scahng theory, 

(e(x)e(0))Gau..ian OC j 1 , (2.9) 

it is possible to work out the following relation between ^ and the thermal critical exponent 
v (we have d = 2): 

^ = 7r(2i/ - 1) . (2.10) 

It yields 

Ml TT , , , , 

— = 2 sin- 2zy - 1 . (2.11) 

Such a result is very important because it gives an exact prediction for the ratio (T2 (T)/o"(T) 
near the deconfining point at Tc as a function of the critical exponent i'. 

Large distance behaviour of correlators: The previous analysis of the mass spectrum 
allows to compute the leading behaviour of the correlators (era) and {VV) at large distance 
by means of their spectral expansion over form factors (the interested reader can refer 
to for the details). On general grounds, the leading behaviour at large distance is 
expected to obey an exponential decay involving the mass of the lightest state allowed by 
symmetry and locality. 

The analysis of the previous section immediately allows to write down the leading term 
for (era) and (VV) correlators in the high-T phase of the theory, up to a proportionality 
constant 

{a{x)a{0)) ~ Ko{M\x\) , |x| ^ oo ; 
{r{x)rm - Ko{Mi\x\) , \x\^oo, 

where Kq denotes the modified Bessel function of order zero, and M, Mi are the masses 
of the soliton and the first breather respectively. It is interesting to notice that the spectral 
expansion gives the exact asymptotic form of the correlator, and not a generic exponential 
decay. 

An important consequence of this observation is that in the regime in which the 
Svetitsky-Yaffe correspondence holds we expect the same large distance behaviour for the 
effective string correction (i.e. the term ^ log R) for both the fundamental and the excited 
string. 

Let us explain this point in more detail. The effective string correction in the case 
of the cylindric geometry of the Polyakov loop correlators has two very different regimes: 
for distances R between the Polyakov loops smaller than L/2 (where L denotes the length 
of the lattice in the compactified time direction) the correction is the usual Liischer term 
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proportional to \/R. On the contrary for R > L/2, it is given by an universal logarithmic 
correction: ilog(^) (see for instance eq. (10) of ref. p8|). 

This is the regime (small L, i.e. high T) in which we may expect the dimensional 
reduction picture to hold, and, according to the identification between 2D spin model and 
3D gauge theory observables discussed above, this implies that the prefactor in front of the 
exponential decay of the two point correlators must be 1 /\fR. This is exactly the prefactor 
of the Kq function and this coincidence represents a non trivial test of the reliability of the 
dimensional reduction program (see the discussion in Sect. 2.2 of [p9| ). 

What is remarkable in the result of Eq. ( p. 12 ) is that this same large distance effective 
string correction holds unchanged both for the fundamental and for the excited string. This 
represents a strong constraint for any consistent effective string model for excited k-strings 
and is one of the exact predictions on the k-string behaviour that we can extract from our 
dimensional reduction analysis. 



Let us summarise the results we discussed in this section in the perspective of applying 
them to the evaluation of the ratio of string tensions near the critical line: 

1. The exact large distance asymptotic behaviour of the correlators of a and V can 
be considered as a reliable tool to extract the lightest mass which governs their 
exponential decay. One can separately compute the Polyakov-Polyakov correlators in 
the representations / and f ® f with a Monte Carlo simulation and then fit the data 
in order to extract M and Mi . Their ratio Mi /M is an estimate of the ratio of the 
string tensions (T2/ cr near the deconfining transition, and can be directly compared 
with its estimate at zero temperature in order to confirm or reject our conjecture. 

2. An independent way to compute the ratio Mi/M is to explicitly use the mass formula 
as a function of the thermal exponent v of the gauge theory which, according to the 
SY conjecture, coincides with that of the corresponding 2D model. In particular one 
can study the finite size behaviour of the plaquette operator (or the susceptibility) 
and extract the corresponding value of v. Then by plugging it in the mass formula 
( |2.11| ) one gets another independent estimate of Mi/M, and again it can be compared 
to the corresponding estimate oia^jo at zero temperature. This is also a direct check 
that the ratio of string tensions follows the proposed analytic formula. 

3. In the large distance regime L < 2R both the fundamental and the excited string 
should be affected by the same effective string correction: ^ log R. 
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2.2 Baryon vertices and mass spectrum 

The general principle invoked in |17] to derive Eq. (1.1) is simply that in a d-dimensional 



gauge theory whatever correlation function made with Polyakov loops in the fundamental 
representation should be described, at sufficiently low temperature and in the IR limit, by 
a two-dimensional conformal field theory with central charge c = d — 2. 

A simple consequence of this general principle is that, as long as the temperature is far 
from the critical one, the shape of the world-sheet spanned by the baryon vertices should be 
temperature independent; this ensures that the baryon static potential has the expected 
asymptotic form |17]. Depending on the location of external sources some fundamental 
strings contributing to the baryon vertex may coalesce, giving raise to the k-string forma- 
tion. As noticed in |17], the balance of the string tensions for a given vertex gives the 
following expression for the angles at the center of the junction of three arbitrary k-strings 

<j%T) + al{T) - a^^{T) ^ ^ 

cos 'di = — — - — - — , and cyclic permutations of the indices. (2.13) 

2aj{T)ak{T) 

The rigidity of the geometry of the vertex is then ensured by requiring that such angles 
are kept fixed when the temperature varies. As a consequence, all the string tension ratios 
are constant up to a given order in T, namely as far as the effective string picture is valid. 
In other words, the previous geometrical construction is likely to break down when the 
system approaches the deconfining temperature, as the string begins to fluctuate wildly. 

We can summarise the above consideration by saying that in the low temperature 
region the trajectory described in the phase space by the gauge system while varying the 
temperature T is a line of constant physics, i.e. a2{T) / a{T) is constant. 

A similar picture emerges when studying the gauge system near the deconfining transi- 
tion. In the framework of the Svetitsky-Yaffe conjecture the second order phase transition 
of the gauge system is described by the critical behaviour of a certain 2D spin model. 
Then, the off-critical scaling region of the latter can be described by a suitable (1-|-1)D 
QFT ^. Focusing on its scattering properties, it is possible to argue that the S-matrix is 
characterised by its analytic properties |31]. In particular the two-particle elastic scatter- 



ing matrix is a multivalued function of the Mandelstam variable s [6 is the rapidity which 
parametrise momentum and energy) 

s = + rn^ + 2mim2 cosh(0i — 82) , (2.14) 

with a Riemann sheet possessing three branch points (it is otherwise meromorphic) , called 
physical sheet. The branch points are at (mi -|- 771-2)^, (mi — 1712)'^ and 00, and the cuts 



■^For a recent application of this approach to 3D Z3 gauge theory see 
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are located on the real line avoiding the interval [(mi — 1712)^, {mi + 1712)'^], which is the 
interval where bound state poles can appear (physical strip). Their masses are eventually 
given by 



mi. 



ml + 777-2 + 2mim2 cos u\2 , triangle of masses , 



(2.15) 



where 6 = i is the purely imaginary value of the rapidity corresponding to the creation 
of the particle m^,. 

Let us stress that this is merely a consequence of the kinematics, and it is a generic 
situation for any (1+1)D QFT. It is nice to see that we are left with the very same structure 
of angles as Eq. ( 2.13| ), which is translated in the so-called "triangle of masses" for the 
bound state particle m^,. Interestingly, since such masses play the role of string tensions, 
we see that they obey the same relation both near T = and near T = Tc, see Fig. |2[ 

Coming back to the case analysed in this paper, we can add a further important element 
to this picture. As explained in the previous Sections, the 3D Z4 lattice gauge theory is 
mapped via Svetitsky-Yaffe to the 2D Ashkin- Teller model near the self-dual critical line, 
which is in turn described by the Sine-Gordon field theory. 

One of the main consequences of the integrability of the latter is the exact knowledge 
of the mass spectrum. In the present case the process of coalescence of two fundamental 
strings into a 2-string corresponds to the scattering of a soliton/anti-soliton pair creating 
the bound state Bi. 



Baryon vertices 




Bound states in integrable QFTs 

mi 



9 = iu 



m2 




"Balance of string tensions" 

, ^ aKT)+'^'P)-al{T) 
cos Ufc — 2cr,(T)cr,(T) 



"Mass triangle" 
m^ = m\ + m\ + 2mim2 cos 



Figure 2: Pictorial representation of the relation between baryon vertices and mass spectrum. 



For such a process we know that = n — ^ which, once inserted in (^.15 ), gives 

M? = 2M2(1 -cosf) ^ — = 2sin-, (2.16) 
which is nothing but the mass formula used in the previous Section. 
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The crucial point is now to point out that the balance of string tensions near T = 
and the mass triangle near Tc have another common feature. In the former the scaling 
properties of the baryon potential require that all the angles 6k should not depend on 
T, therefore the string tension ratios (Jk/cr stay constant as T varies and define a line of 
constant physics starting at T = 0. In the latter the angles involved depend only on the 
marginal coupling ^ hence a variation of the relevant coupling r of the 2D model ( |2.2| ) 
generates a line of constant physics starting at r = 0. 

In the gauge/CFT map established by the SY conjecture, in order to avoid a mismatch 
between the RG trajectories of the 3D gauge system and the corresponding 2D model, 
relevant perturbations of the CFT should correspond to relevant couplings of the critical 
gauge system. On the gauge side we pass from the T ~ region to the T Tc by keeping 
constant the gauge couplings (hence also the lattice spacing) and varying simply the size of 
the imaginary time direction. Thus the only physical parameter which is varied in passing 
from low temperature to Tc is the reduced temperature t = '^^^ of the gauge system, 
hence near r ~ we have r = T(t), while ^ is kept constant. 

Summing up, the variation of temperature of the gauge system defines a line of constant 
physics near T = and a similar line near T = T^- The numerical work in the next two 
Sections will show that these two lines are in fact a single one which goes through the 
whole confining phase. 



3. Monte Carlo setting and procedure 
3.1 Mass ratio by correlators 



As introduced in Subsection 2J , we can determine the ratio M\ /M using the large distance 
asymptotic behaviour of correlators; actually, exploiting the Svetitsky-Yaffe conjecture, we 
measured the Polyakov-Polyakov correlators G^{R) of the (2+l)D Z4 gauge theory: 

G^{R) = (P^(O)Pi(i?)) . (3.1) 

In Section we have explained we can study this theory by means of simulations on the 
dual 3D AT model and in [^2| the measurement of Polyakov-Polyakov correlators in both 
the fundamental and double fundamental representations, G{R)f and G{R)ff, is described 
in detail. 

Note that here, since we are studying the theory near the critical line, the periodic 
boundary conditions play an important role; in this case there are 16 topologically different 
surfaces bounded by the two Polyakov lines (4 for each Ising variable, as discussed in p^]), 
therefore the correlator is the sum of these contributions. In fact, there are only two 
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important contributions, so we take into account only these two in our simulations (see 
Figures ^ and |^ . 
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Figure 3: Polyakov-Polyakov correlator in the fundamental representation. 
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Figure 4: Polyakov-Polyakov correlator in the double fundamental representation. 

We have taken 10^ measures on the 64? x 7 lattice in the phase space point P; N-r = 7 
is chosen because it is the lowest possible value above the deconfinement transition. Points 
on the plots are obtained by independent simulations, one for each value of R in the 
range [15 -r- 44]. These data are fitted using an expansion of the Ko{mR) Bessel function, 
truncated to first two terms, 

mi? 



G{R) = const X 



\/mR 



1 



1 

8mR 



+ "echo terms". 



(3.2) 
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in a range [Rmin, Rmax], where Rmax = 44; we have verified the results are stable when 
R-min varies in the range [22 33]. Therefore, it is possible to determine the two masses: 

aMff = 0.0698(15) (x^/d.o.L w 1.3) , 
aMf = 0.0433(8) (xVd-o.f. ^ 1.2) , (3.3) 

from which we can determine the ratio: 

o-2(r - Tc)/a{T ~ Tc) = Mff/Mf = 1.612(46) . (3.4) 

This result, obtained near the critical temperature, is compatible with the zero-temperature 
value ( |1.10D , providing a strong evidence for our conjecture. 

3.2 Estimating a2/cr through the thermal exponent v with finite-size scaling 



To use the formula for the mass ratio, Eq. ( 2.11| ), we need a quite precise estimate for the 



thermal critical exponent v in the phase space point P (see Fig. ||). It can be obtained 
by means of a finite-size scaling analysis once the critical temperature for that choice of 
couplings has been identified. 

The problem is, the system at the coupling P turns out to be critical for a temperature 
Tc such that 6 < < 7, hence, having to work with integer inverse temperatures, it is 
not possible to avoid some approximate method. The idea is then the following: choose a 
direction in the (a, f3) phase space which crosses the critical line, and by moving from P (in 
opposite directions) along it find two new points, Pf and Pq, at which the system is critical 
for temperatures T = 1/7 and T = 1/6 respectively. There, perform an estimate for v 
with a standard finite-size scaling approach. Then, with a linear interpolation, construct 
the corresponding quantity for the original P (see Fig. ^). 

This first-order approximation, however, is motivated only if the points are close 
enough that the variation of v is quite small; since the exact shape, in the phase space, of 
the trajectories of the RG (that is, the set of points which are mapped to the same point 
on the Gaussian model) is unknown, the best guess is to move from P along a direction 
which is perpendicular to the zero-temperature critical line connecting the two decoupled 
Ising and the 4-state Potts systems. 

It has been shown long ago that in the SY context the plaquette operator is mapped 
into a combination of the unity and the energy operator of the corresponding CFT p^ . 
Therefore, once the system is made critical, one could extract the thermal exponent v from 
the finite-size scaling behaviour of the plaquette operator or some related observable that 
we denote with ( □ )l) where L is the spacial size of the lattice. 
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Figure 5: In the three-dimensional AT system (on the left), the point P belongs to a whole 
line of points which exhibit the same critical behaviour, which is a RG trajectory in phase space 
(represented as a solid line passing through P). The intersections of this line with the finite- 



Relying on a linear approximation, however, we located the points P6, P7 as described in the text, 
in order to identify, by interpolation, the point on the dimensionally-reduced AT model (on the 
right) to which P is mapped. 

In order to exploit the computational advantages of the dual transcription of the gauge 
model, it is convenient to evaluate directly the internal energy of the 3D AT model defined 
in (|1.6| ), namely 



We expect the foUovifing finite-size critical behaviour as a function of the spacial side L of 
the system: 



temperature critical lines for ^ = 6, 7 are labelled Q6 and Q7, but their exact position is unknown. 





(3.5) 




(3.6) 



If one, instead, looks at the corresponding (density of) susceptivity: 




(3.7) 



the power-law to compare with has the form: 





with the advantage that no constant additive terms are present, which could largely spoil 
the stability of the numerical results. 



To locate the critical points Pq and Pj, we proceeded in the phase space in a dichotomic 
way along the above-mentioned line from P and looked for peaks in the plaquette suscep- 
tivity on a fixed spatial size system. There is an intrinsic uncertainty on the exact critical 
values for a, (3, but in principle it can be indefinitely shrunk by considering larger and 
larger lattices. 

To perform the simulations, we used a cluster-based nonlocal update algorithm, an 
adaptation of the Swendsen-Wang prescription, based on alternating global updates on the 
two Ising subsystems in which the other variables play the role of a frozen background 
field. The algorithm is described in more detail in [p^]. 

We used L = 200 finite-temperature lattices to find the couplings corresponding to Pg 
and P7, and for each sampled value of the couplings we took at least 0(10^) measurements. 
By locating the peak in the plaquette susceptibility (see Fig. we could identify the two 
points with a certain degree of accuracy as Pe (0.0500965 ± 0.0000063, 0.207235 ± 0.000015) 
and P7(0.0497859 ± 0.0000077,0.206478 ± 0.000019). 
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Figure 6: Behaviour of the plaquette susceptibility in the AT model at T = 1/7 and spatial size 
L = 200. The clear peak allowed a precise estimate of the critical point. 



Then, on exactly critical systems at Pq and P7, we took 0(10^) measurements of the 
plaquette at 26 values of spatial side L, ranging from L = 10 to L = 165. Not surprisingly, 
from the bare plaquette data the signal of the power-law behaviour was very noisy due to 
the presence of a constant background as another fit parameter, so we switched to using 



the prediction (3^) for the susceptibility. 

The data fitted very well to the expectation from L = 70 already, so we could extract 
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two values of the critical index v (see Fig. |^: 

vt^^Iq = 0.8004(19) [22] ; 

i^T^^/j = 0.7942(18) [38] , (3.9) 

in which the first uncertainty refers to the statistical fluctuations while the second is an 
estimate of the systematic error in the measurement. 
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Figure 7: Finite-size scaling of the susceptibility of Eq. i.i in the AT model at T = 1/7 as a 
function of the system spacial size. The data allowed a precise estimate of the critical index v. 

By linear interpolation along the couplings, the value of u and the (coupling-dependent) 
critical temperature Tc was calculated for the very point P. It was found that Tc{P) — 
0.1502 ~ 1/(6.655), which (since the value of a is well known for P) gives the universal 
ratio 

T. 

= 1.0393(12) . (3.10) 



From the interpolation, we have i^(-P) = 0.7984(19) [27]. By plugging it into the formula 
for the mass ratio ( |2.11| ), we obtain the following result: 

Ml 



M 



(P) = 1.6124(71)[102] 



(3.11) 



which is compatible with the less accurate estimate coming from the quantities in [|17| and 
thus well supports our conjecture. 



4. Conclusions 

In this paper we studied the ratio of the string tensions a2{T)/a{T) near the deconfining 
point Tc of a 3D Z4 gauge model and compared the result with a general formula which is 
expected to be true near T = for a generic gauge theory in three or four dimensions. 
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In this particular case we have combined numerical experiments with known exact 
results of an integrable 2D quantum field theory that belongs, according to the Svetitsky- 
Yaffe conjecture, to the same universality class of the critical gauge system. 

An interesting property of the integrable model is that the mass ratio of the two 
physical states of the theory, which should equate the string tensions ratio near Tc, can be 
expressed as a simple function of the thermal exponent v (see Eq. ( 2.11| )). Therefore we 
used two different methods to evaluate such a ratio: either a direct evaluation of the string 
tensions through Polyakov loop correlators near (see Eq. ( |3.4| )) or through a measure of 
V (see Eq. ( |3.11 )). Both the estimates give compatible results which nicely agree with the 
ratio (72 /o" evaluated at T = (see Eq. ( |1.1[)| )); for a schematic summary, see Table |l|. We 
then conclude that, at least in this model, the k-string tensions ratios do not depend on T . 



^2(T)/a(r) 


Temperature Method 


1.610(13) 
1.612(46) 
1.6124(71)[102] 


r~o + o(r^) Eq. (i.lop 

T < Eq. (P) 
T < Eq. (3.11D 



Table 1: Numerical results. 



Even if in Section 2.2 we gave a general RG argument to support this assumption in a 
wider context, we do not dare to extend such a conjecture to a general gauge system, one 
reason being that if the deconfinement transition is first order, as is the case in most gauge 
theories, we do not know a sound argument to support it. 
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